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Abstract 

We have studied the relevance of spin-orbit coupling to the dispersion relation of the Larmor 
resonance observed in inelastic light scattering and electron-spin resonance experiments on GaAs 
quantum wells. We show that the spin-orbit interaction, here described by a sum of Dresselhaus 
and Bychkov-Rashba terms, couples Zeeman and spin-density excitations. We have evaluated its 
contribution to the spin splitting as a function of the magnetic field B, and have found that in the 
small B limit, the spin-orbit interaction does not contribute to the spin splitting, whereas at high 
magnetic fields it yields a B independent contribution to the spin splitting given by 2{X 2 R — \ 2 D ), 
with Xr^d being the intensity of the Bychkov-Rashba and Dresselhaus spin-orbit terms. 

PACS numbers: 73.21.-b, 73.21.Fg 
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I. INTRODUCTION 




The study of spin-orbit (SO) effects in semiconductor nanostructures has been the object 
and theoretical investigations in the last few years, see e.g. Refs. 
I2I13B14I and Refs. therein. It links the spin and the charge dynamics, 
hence opening the possibility of spin control by means of electric fields . 15lie 

It has been recently shown^ that the SO interaction affects the optical properties of 
GaAs quantum wells by inducing a coupling between charge density and spin density ex- 
citations in the long wavelength limit. We extend here this study to the influence on the 
Larmor resonance of the combined effect of both Dresselhaus^ and Bychkov-Rashb a 19 ^ 20 SO 
interactions, and use our results to discuss some features of the spin modes disclosed by 
inelastic light scatterin g 21 ' 22 and electron-spin resonance experiments . 2 ? 124 

Our approach is based on the solution of the equation of motion up to second order 
in the SO intensity parameters.— This method has been also used to derive the Kohn 
theorem,— and goes as follows. We write the Schrodinger equation for a iV-particle system 
as H\n) = E n \n), with |0) and E being the ground state (gs) and gs energy, respectively. 
If one can find an operator 0+ such that |n) = O+|0), O n \0) = 0, it is possible to cast the 
Schrodinger equation into an operator equation -the equation of motion- [H,0+] = ^ n O^, 
where uj n = E n — Eq is the excitation energy of the state |n). The solutions of this equation 
are used to find the excitation energies of the system as well as its excited states in terms 
of their creation operators. 

This work is organized as follows. In Sec. II we apply the equation of motion method 
to the Larmor mode in the presence of a SO coupling. The results are used in Sec. Ill to 
discuss the spin modes in quantum wells, and are compared with the experimental results 



of Refs. 
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II. THE EQUATION OF MOTION APPROACH AND THE LARMOR MODE 

The operators describing the SO Rashba and Dresselhaus interactions are respectively 
given by 

A N 

11 " (!) 

h 3=1 
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and 

A N 

H D = fY,[ p ^- p y°yh > ( 2 ) 
i=i 

where the ex's are the Pauli matrices and P = — ihV+-A represents the canonical momentum 
in terms of the vector potential A which in the following we write in the Landau gauge, 
A= B(Q,x,0), with B=Vx A = Bz. 

In the effective mass, dielectric constant approximation, the quantum well Hamiltonian 
H can be quite generally written as H = Hks + V res , where Hks is the Kohn-Sham (KS) 
one-body Hamiltonian consisting of the kinetic, Rashba, Dresselhaus, exchange-correlation 
KS potential and Zeeman terms, and V res is the residual Coulomb interaction. The KS 
Hamiltonian reads 



N 



Hks = 



P+P- + P-P+ \ R ,^ „_ , A; 



-(P+a_ - p-a+) + -^{P + a + + P-a-) 



Am 2ih y ' 2Ti 

(3) 



+W xc (n, C, V)a z + -g*fi B Ba z 



^ 
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where m = m*m e is the effective electron mass in units of the bare electron mass m e , P ± = 
P x ± iP y , and cr± = a x ± io y . Although other approaches may be also considered, we have 
considered the exchange-correlation potential W xc (n,£,V) in the local-spin current density 
approximation (LSCDA) . 26 i 27 It depends on the density n, magnetization C, = — n^, and 
local vorticity V, and is evaluated from the exchange-correlation energy per electron S xc as 
W xc = d(nS xc ) I 'd£. The last term in Eq. Q is the Zeeman energy, where //# — ^e/ (2m e c) is 
the Bohr magneton, and g* is the effective gyromagnetic factor. For bulk GaAs, g* = —0.44, 
m* = 0.067, and the dielectric constant is e = 12.4. To simplify the expressions, in the 
following we shall use effective atomic units h = e 2 /e — m — 1. 

In the following, the residual Coulomb interaction will be treated in the adiabatic time- 
dependent LSCDA (TDLSCDA)j2Z. We are going to see that, in the absence of SO coupling, 
not only the exact Hamiltonian, but also the one in which the residual interaction is treated 
in the TDLSCDA fulfill the equation 

[H,S T ] = ±u L S T , (4) 

where S T = an d uj l = \g*fi B B\. Thus, if |0) is the gs of the system, the states 

S T \0) are eigenstates of H with excitation energies ±o>£. Note that a negative g* implies 



that the spin-up states are lower in energy than the spin-down ones, and that the actual 
physical solution of Eq. (£Q) is that corresponding to the S- operator. This is the physical 
contents of the Larmor theorem. Note also that in the absence of spin-orbit coupling, 
[H, Pj~] = u c J2f Pj~, where cu c = eB / (mc) is the cyclotron frequency. This is the Kohn 
theorem, which also holds in the adiabatic time- dependent local spin density approximation 
(TDLSDA) and in the TDLSCDA, and can be generalized to the case of quantum wires and 
dots parabolically confined. 
Since 

N 

[H, S-] = lu l S- + 4^[A d P + ( x z + i\ R P-a z } 3 , (5) 
i=i 

the spin-orbit terms in Eq. (J3j) mix the transverse spin excitations induced by the operator 
S- with the spin-density excitations induced by J2jLi Pj ^ii an d thus Larmor's theorem is 
not fulfilled. In the following, we use the equation of motion approach to find the eigenvalues 
and eigenstates of the KS Hamiltonian Hks Eq. (jBJ) which arise from the SO mixing, and 
will evaluate the spin wave dispersion relation w(q) by taking into account the effect of the 
residual interaction. This is done by firstly solving the equation of motion 

[H KS ,0 + ] =uoO + , (6) 

and then calculate the transverse response Xtiq, per unit surface A in the TDLSCDA: 



where F xc = W xc /Ci an d X^iQ^) is the KS transverse response per unit surface^ 27 The 
poles of Xt(q, yield ui{q). The transverse spin response without inclusion of spin-orbit 
coupling has been studied in the past in the RPA 28 and time- dependent Hartree-Fock 29 
appr oximat ions . 

Up to second order in A^^, Eq. (jUJ) is straighforwardly solved by the operator + = 
Y.j I =i[ a<J - + bP + a z + cP~a z + da + ]j. To do so, one has to use the commutators [a + , aJ\ = Aa z 
and [P~ , P + ] = 2u c . This yields a homogeneous system of linear equations for the coefficients 
a, b, c and d from which the energies to are obtained by solving the secular equation valid 
up to \ 2 R D order 

(to 2 - u 2 L ){u 2 - cu 2 ) -4lu c (\ 2 d + \ 2 r )lu 2 -4lu 2 Cu l (\ 2 d - X 2 R ) = 0, (8) 
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where ujl = \g* I^bB + 2W XC \. This quartic equation can be exactly solved, yielding the 
excitation energies (only positive solutions are physical). For each of them, the homoge- 
nous linear system, supplemented with the normalization condition (0|[(O + )^, O + ]|0) = 1, 
determines the coefficients a, b, c and d. We have found it more convenient to discuss the 
solutions of the above equation in the limits of small and large magnetic fields, which are 
more transparent and easier to compare with available experimental data. In the small B 
limit (ujl, uj c <C Xr, Ad) we obtain a unique solution 

00 = 2 y^ c (A| + XI) . (9) 

In the large B limit we obtain 

cu(S T ) =±(u L + 2X1-^- - 2A D — , (10) 

V OJc + ^L UJ C -UJ L J 

which are mainly excited by the operators S*_ and S + , and 

wCP^) =±U- 2X1^^- + 2X1—^-) , (11) 

V U) c + UJ L UJ C -UJ L J 

which are mainly excited by the operators X/^Lit-P 4 'o ' z ]j and J2jLi[P~^z]j- By mainly we 
mean that the coefficient of the corresponding operator entering the definition of is 0(A^ d), 
whereas all the others are 0(X 2 RD ). Note that if A^d = 0, the two physical modes in the 
preceding equations are uncoupled. 

Equation Q shows that, at B ~ 0, to order X 2 RD there is no spin splitting due to the 
SO coupling. Indeed, when B — > not only lvl and uj c vary linearly with B, but also W xc 
does, implying that the solution of Eq. (JHJ) goes to zero in this limit. Earlier electron-spin 
resonance measurements on GaAs quantum wells^ seemed to indicate that a finite spin 
splitting was present in the B = limit. However, subsequent experiments carried out by 
the same group^ covering a broader B range point out that the spin splitting of a Landau 
level is an exact quadratic function of B, and that its extrapolation to B = leads to 
a vanishing spin splitting. Our result, which is not changed by the effect of the residual 
interaction, is thus in full agreement with the experimental findings of Dobers et al.— 

We have checked that, at low B fields, the dominant component of + corre- 
sponding to the energy Eq. (JH1) is the spin-flip operator Y J 1 j=i\ cr -}j- ^ n this limit, 
the Dresselhaus and Rashba SO interactions act "in phase", whereas at high 
B they partially compensate each other [compare the energy given in Eq. d^I) 
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with these in Eqs. (|10|) and This arises from the structure of the secular 

Eq. (JHj), where in the low B limit, the second term dominates over the third 
one, whereas in the high B limit both terms are equally important, yielding 
the solutions shown in Eqs. (|10|) and We have not been able to find a 

deeper explanation for this different behavior at low and high magnetic fields. 
It is worth to mention that the independent particle Hamiltonian can be exactly 
solved when only the Rashba or Dresselhaus SO terms are included. 30 The merit 
of Eqs. (pi llO|) is that they are exact to the relevant \ 2 R D order when both SO 
couplings are simultaneously taken into account. 

The excitation energy uj(SJ) is the independent particle (KS) value for the spin splitting 
and violates Larmor's theorem even if the SO coupling is neglected. On the contrary, when 
the residual interaction is properly taken into account, the theorem Eq. (J3J) is recovered. In 
the following, we will concentrate on the large B limit. It is then possible to derive the spin 
wave dispersion relation, including spin-orbit effects, by solving the equation 

l-2F xcX ? S (q,uj) = (12) 

that gives the poles of the transverse response function Eq. (J7J). To do so, we write the 
transverse spin response as 27 

AXt (9.w) = J l — , 13 

uj — u(S-) u> + u>(b + ) 

where \u>(S T )) = O + [u(S T )]\0) , and the corresponding energies are given by Eqs. (fTUj) and 
(111)1 ■ The calculation of the matrix elements in Eq. (|T3*|) must be done with care since |0) 
and \u)(S±)) are not eigenstates of S z because of the spin-orbit coupling. Neglecting terms 
in X 2 rd /uj c or smaller, one gets 3 ^ 

^')=«3j' (14) 

where 

1 N 

n?) = Tf(0|E^ rj |0) (15) 

iV 3=1 

is the gs elastic form factor. Note that F(0) = 1 and that F(q) goes to zero when q — > oo. 
From Eq. ()12)) one finally obtains 
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u = \g*ti B B\ + 2X 2 R ^^- - 2X1^^- - 2W XC (l - |F(g)| 2 ) (16) 

This is the main result of our work, together with the lack of SO splitting we have found 
in the small B limit. In the limit q — > oo, Eq. (fT7)J) yields the independent particle spin 
splitting Eq. (|TU|) . which crucially depends on the actual value of W xc entering the definition 
of Ul- In the q = limit, neglecting terms of order lj l /lj c , Eq. (JT5j) reduces to^i 

u = \g*Li B B\ + 2(\l-\ 2 D ). (17) 

This expression shows that, neglecting the SO coupling, Larmor's theorem is fulfilled in the 
adiabatic TDLSCDA (it can be shown that the same holds in the adiabatic TDLSDA), and 
that at high magnetic fields, the SO interaction yields a B independent contribution to the 
spin splitting. Taking e.g. m\ 2 R /Ti 2 = 27 /zeV, m\ 2 D /Ti 2 = 6 /xeV, which have been recently 
used to reproduce the spin splitting in quantum dots^ 4 - and the splitting of the cyclotron 
resonance in quantum wellSj 17 we get 2m(\ R — \ 2 D )/h 2 ~ 40 /xeV. This is definitely a small 
amount, but it may have an influence on the fine analysis of some experimental results (note 
the vertical scale in Fig. QJ. 



III. COMPARISON WITH EXPERIMENTS AND DISCUSSION 

Using inelastic light scattering, Davies et al.— and Kang et al.— have measured charge 
and spin density excitations in 2D electrons systems confined in GaAs quantum wells at 
high B. In the following we only discuss the results of Ref. |2l| because the information 
presented in Fig. 2 of this reference is especially well suited for the purpose of our work. 
These results are represented in Fig. ^ The data labeled S correspond to wave-vector 
allowed scattering from the q = Larmor mode; indeed, the maximum in-plane q allowed 
by the experimental geometry is small, q ma x = 6 x 10 4 cm -1 , so that the S mode of energy es 
should correspond to the spin splitting energy Eq. (fTTj) . The data labeled SW is attributed 
to disorder- activated scattering, and would correspond to q ^ excitations of energy egiy.— 
The difference between e$w an d £s is attributed to the exchange enhancement of e$w above 
the Zeeman energy.— The dashed straight line represents the Larmor energy taking for g* 
the bulk value, \g*\ = 0.44. This overestimates es, especially at high B. 
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Any sensible comparison with these results must take into account the B dependence of 
g*. This dependence has been clearly established in magnetoresistivity experiments,— taking 
advantage of the fact that the electric-spin resonance affects the magnetoresistivity of the 
2D electron gas, and this can be used to determine g*{B). These experiments probe the 
one-electron energy levels and are not influenced by many-electron interactions, contrarily to 
magnetoquantum oscillations, which are strongly influenced by many-electron interactions. 



The spin splitting obtained in Ref. 



24 is represented in Fig. as dots and crosses, which 



correspond to two different samples. In the lowest Landau level, which is the common 
physical situation for all data represented in Fig. ^ these authors have fitted g*{B) as 
\g*(B)\ = #q + rB/2. The value of the parameters r and g$ turns out to sensibly depend 
on the experimental sample, and the possibility of a SO shift at high B values could not 
be considered there. Moreover, the A R y n values are rather poorly known and dependent on, 
e.g., the thickness of the experimental sample. 



We have thus renounced to use the laws g*(B) obtained in Ref. |24| in conjunction with 
Eq. (|17jl. to establish a clear evidence of spin-orbit effects on the es energy obtained from 
resonant inelastic light scattering experiments, and have satisfied ourselves with the more 
limited scope of using Eq. (jl7j) as a three-parameter law to fit es as well as the spin splittings 
of Ref. 13, with the aim of seeing whether a reasonable value for these parameters can be 
extracted. 

The solid straight lines in Fig. represent the result of such linear fits, whose parameters 
are collected in Table HI In the case of inelastic light scattering, the neglect of the SO term 
in Eq. (fTTj) yields an unrealistic g$ = 0.49, as this value should be smaller than that of bulk 
GaAs due to the penetration of the electron wave functions into the Al^Gai-^As barriers. 
The dispersion of the electron-spin resonance datapoints^ seems to be smaller, and the 
analysis of the high B data might be used to ascertain which SO mechanism is dominating 
in a given sample. This could be an alternative or complementary method to the recently 
proposed^ of using the anisotropy of the spin splitting in single-electron resonant tunneling 
spectroscopy in lateral quantum dots submitted to perpendicular or parallel magnetic fields. 
The analysis of samples 1 and 2 would indicate that in the former, the Dresselhaus SO is the 
dominating mechanism, whereas in the later it is the Bychkov-Rashba one. We want to 
stress that we have extracted the experimental data from a careful digitalization 
of the original figures. Due to the smallness of the effects we are discussing, we 
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cannot discard that this procedure may have had some effect on the value of the 
parameters determined from the fit, and our analysis should be considered as 
qualitative to some extent. However, we find it encouraging that the parameters 
obtained from the fits are meaningful, and within the range of values found in 
other works. 1 &£L 

We finally discuss briefly the q ^ SW mode. From Eqs. ()16j) and ([17)1 . we have that 
£sw — £s = —2W XC (1 — \F(q)\ 2 )- At high q, this difference is sensibly determined by W xc . 
In our calculation, as well as in time-dependent Hartree-Fock 2 ^ and exact diagonalization^ 2 . 
calculations, we have found values of W xc of the order of —2 meV. Hence, —2W XC is about 
a factor 40 larger than the measured esw ~ e s, which has been obtained in the q — > 
limit where short range correlations are very important in determining the actual value of 
1 — |F(g)| 2 . This can be seen by assuming for F(q) the independent particle value. Using 
a Slater determinant made of Fock-Darwin single particle wave functions to describe the 
gs of the system at B ^ 0, one finds F(qi) = e _?2£2 / 4 , where £ is the magnetic length 
i = (hc/eBf/ 2 . In the small q limit, e sw -e s ^ -W xc q 2 £ 2 . For B = 10 T, q max i ~ 0.05 and 
tsw — ~ 0.01 meV, which is about one tenth of the experimental result as shown in Fig. 
^ Light scattering experiments at small q and high B are thus very sensitive to correlation 
effects in the elastic form factor, which is the key quantity to reproduce the experimental 
findings. 
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TABLE I: Parameters of the least-square linear fit to the cited experimental data. 





h 2 m(X 2 R - A 2 n ) 

V it L) / 

(/ieV) 


* 


r 

(T- 1 ) 


e s (Ref. 21) 


8.6 


0.43 


-7.3 x 10" 3 


samplel (Ref. 24) 


7.7 


0.38 


-8.1 x 10~ 3 


sample2 (Ref. 24) 


-1.3 


0.41 


-1.3 x 10" 2 
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FIG. 1: Experimental spin splittings as a function of B. Squares and diamonds, from Ref. |21 
dots (sample 1) and crosses (sample 2), from Ref. I24. The dashed line corresponds to ^b9*B with 
\g* \ = 0.44, and the solid lines are fits using the law Eq. (|17|). 
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